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Abstract 

It was recently suggested by Blythe and Evans that a properly denned steady 
state normalisation factor can be seen as a partition function of a fictitious statistical 
ensemble in which the transition rates of the stochastic process play the role of fu- 
gacities. In analogy with the Lee- Yang description of phase transition of equilibrium 
systems, they studied the zeroes in the complex plane of the normalisation factor in 
order to find phase transitions in nonequilibrium steady states. We show that like 
for equilibrium systems, the "densities" associated to the rates are non-decreasing 
functions of the rates and therefore one can obtain the location and nature of phase 
transitions directly from the analytical properties of the "densities" . We illustrate 
this phenomenon for the asymmetric exclusion process. We actually show that its 
normalisation factor coincides with an equilibrium partition function of a walk model 
in which the "densities" have a simple physical interpretation. 

1 Introduction 

The extension of concepts used in equilibrium statistical mechanics, like the free energy, 
to nonequilibrium steady states has a long history ^H]. That a simple extension is not 
possible can be seen in [2] where it was shown that in certain cases the free energy func- 
tional is not a convex function of the density. On the other hand Arndt [1J has shown in an 
example that applying the Lee- Yang description using the zeros of an ad hoc definition of a 
grand-canonical partition function gives the correct phase transition. Further applications 
of this idea can be found in 00131 • In a very interesting new development Blythe and 
Evans [H] considered the normalisation of the stationary state of several stochastic systems 
as a function of the transition rates and applied the Lee- Yang approach in the same way 
as one would for an equilibrium partition function. While a normalisation may seem to 
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be defined ambiguously, it was noted in [HlllUj that in a formal way, a unique definition of 
this normalisation can be made using the matrix-tree theorem, which has a long history in 
graph theory going back to Sylvester see also |13U12j and Section |2 This connection 
explicitly relates the normalisation of a stationary state to the combinatorial problem of 
counting weighted spanning trees on graphs, which implies a direct interpretation of the 
normalisation as a statistical mechanical partition sum. Can we can learn anything about 
the steady state phase diagram from the spanning trees? There are a few problems with 
that. For example, it is generically unclear what the correspondence is between original 
dynamical quantities and those that are natural for describing the spanning trees. Sec- 
ondly, the normalisation as defined above may not be "minimial" in the sense that it may 
contain an overall nontrivial polynomial factor which is common to each of the stationary 
state weights. Such a polynomial factor cannot contribute to the phase behaviour of the 
stationary state. We will call the normalisation as defined via the matrix-tree theorem 
but with common factors removed the reduced normalisation. 

The purpose of this paper is to try to bring a better understanding of the Blythe-Evans 
approach which is summarized in Sectional In this section we also show that to each tran- 
sition rate one can formally associate "particle numbers" the same way one relates particle 
numbers to fugacities. Moreover, like in thermodynamics, one can prove that the "particle 
numbers" are non-decreasing functions of the transition rates. This important observation 
allows us to detect the existence of phase transitions from the behavior of the "particle 
numbers" in the space of transition rates. At this point the physical meaning of the "par- 
ticles" is completely obscure (there are as many "particles" as the number of independent 
transition rates minus one). Moreover, the thermodynamic potential defined through the 
(reduced) normalisation factor is not an extensive quantity. The volume, defined by the 
leading asymptotic behaviour of the normalisation, and therefore the definition of the 
"densities" might change in the space of transition rates. This allows for phase transitions 
not encountered in equilibrium models with local interactions. This phenomenon appears 
in the following way. In a certain domain of the fugacities, the "densities" span the entire 
interval between zero and one. This defines a "phase" (inside this domain one can have, 
like in equilibrium, several phase transitions). The boundary of the domain separates it 
from another domain ("phase") where one has to take another definition of the "densi- 
ties" because of the change of the volume. In this second phase the "densities" are not 
necessarily finite. 

The fact that the reduced normalisation factor might have a direct physical interpre- 
tation is known from the raise and peel model This is a one-dimensional stochastic 
model with nonlocal transition rates and its reduced normalisation factor (whose log- 
arithm is proportional to the square of the system size) coincides with the number of 
configurations of the two-dimensional ice model with domain walls boundary condition — 
an equilibrium problem. It is our aim to show that this situation is more general. 

In Section El we define the one-transit walk model (OTW). This model, which is not 
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parity invariant, depends on two parameters which are the Boltzmann weights or fugac- 
itites of contact points. We compute the partition function of this model as well as the 
two densities corresponding to the two fugacities. The two densities have a clear physical 
meaning. The phase diagram of the OTW model is obtained from the expressions for 
these densities. It is the same as that of the totally asymmetric simple exclusion process 
(TASEP) ^JE3E1E3E3 if we replace the two fugacities with the boundary rates of 
TASEP. As we are going to show, the derivation of the phase diagram of the TASEP 
model from the analytic properties of the number of "particles" as a function of fugacities 
will give a better understanding of the nature of the phase transitions. In Section 13.41 
we discuss the microscopic properties of the OTW model making cleare the connection 
with the TASEP. We conclude with a discussion of the phase diagram obtained from the 
"densities" of the partially asymmetric simple exclusion process (PASEP) |2UU21) . In this 
case we have three kind of "particle" numbers: two associated with the boundary rates 
and one associated to the back hopping rate q. At the symmetric point q = 1 a new kind 
of phase transition occurs. 

Our conclusions are presented in Section 0J 



2 The normalisation as a positive polynomial 

Let us start by considering an arbitrary Markov process in continuous time on a state 
space spanned by the states {|a)}" =1 , whose master equation is given by 

^Pt(a) = ]T (r ab P t (b) - r ha P t {a)) . (1) 

The r a h are the transition rates from state \b) to \a) and Pt(a) is the (unnormalized) 
probability to find the system at time t in state \a). Equation Q can be conveniently 
rewritten as 

d n 

-\P t ) = -H\P t ), \P t )=Y,Pt(a)\a), (2) 

a=i 

where H is the matrix with off-diagonal elements H ab = —r ab and whose columns add up 
to zero. One of the main properties of interest of such a Markov process is its long time 
behaviour. In the limit t — > oo the system approaches its stationary state |Poo)> which we 
will assume to exist and for simplicty to be unique, given by 

mPoo) = o. (3) 

The stationary state is thus given by the right eigenvector of the matrix H corresponding 
to its eigenvalue 0. This equation can be solved in the following formal way, see e.g. |22j . 
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Let H (a, b) be the matrix corresponding to H with the ath row and 6th column removed. 
The cofactor X(a, b) is then defined by, 

X(a, b) = (-l) a+b det H(a, b). (4) 

If the eigenvalue is unique, 

= det# = Y,H ab X(a,b) = Y,H ab X(b,b), (5) 

b b 

where we have used X(a,b) = X(b,b) for all a (see Appendix |A")) . We see that the 
eigenvalue equation (j2J) is solved by the cofactors of H, 

H\P)=0, P(b)=X(b,b). (6) 

This solution fixes a particular normalisation of the eigenvector for all system sizes. This 
normalisation is uniquely defined up to an overall rescaling of H, or equivalently a rescaling 
of time (which can vary with the system size). To be able to interpret P(b) as a probability 
distribution, we write 

n n 

Poo(b) = P(b)/Z n , Z n = Y J P{b) = Y, X ^ b )- ( 7 ) 

6=1 6=1 



It can be shown using the matrix-tree theorem |13|I12| that the normalisation Z n of 
a stationary state of any stochastic (Markov) process is always given by a homogeneous 
polynomial in the rates r ab (some of which might be equal or be zero), of degree n — 1 and 
with positive coefficients, i.e. it has the form of a generating function. A simple proof of 
this important statement is for example given in ?22, , which we have included in Appendix 

E3 

We would like to identify the rates as generalized Boltzmann factors or fugacities 
fab = z a b , and Z n ({z a b}) as a generalized partition sum for nonequilibrium systems. Since 
the normalisation Z n is a polynomial in the variables z a \, with positive coefficients, by the 
Cauchy- Schwartz inequality its negative logarithm, 

F n = -logZ n , (8) 

is therefore a convex function in all its arguments z^- In analogy with equilibrium statis- 
tical mechanics, we will associate to each rate r a ^ = z a b a "particle number" iV a b, 

N ab = -z ab ^- (9) 
dz ab 

These numbers are positive and increasing functions of the fugacities for any size of the 
system but they are linearly dependent. One can arbitrarily choose one rate equal to one 
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which fixes the time scale and leaves the remaining rates dimensionless. In this way we 
are left with one fugacity less. In the large n limit, 

N ab = V(n) Pab , (10) 

where V(n) is the volume and p ab are the "densities". One can now use the equilibrium 
approach to the theory of phase transitions and apply it to the densities p ab . A first order 
phase transition, for example, is a location in the transition rates space where one of the 
"densities" has a discontinuity. 

It may happen however, as it will in our example below, that all cofactors X(b,b) 
contain a common nontrival polynomial factor. Such a common factor will cancel out 
in Poa(b) and hence cannot contribute to the nonequilibrium phase behaviour. In © 
however it could give rise to spurious singularities that are not related to the physical 
phase transitions. In the example of Section |21 no such spurious phase transitions appear 
(see Section 13. 3. 2 j) . 

There is another major difference however between equilibrium systems with short- 
range interactions and the present problem: the "particle numbers" are not necessarily 
extensive quantities (see the examples in Sections 13.21 and 13. 5|) . This implies that in the 
parameter space the p ab might diverge and we have to change the definition of the factor 
V(n) in (|1U|) . Actually such a phenomenon is also known in equilibrium problems with 
non-local interactions (see |23j ) in the theory of special surface phase transitions 24 . As 
we are going to show in Section 1331 the analogy goes deeper. 

The philosophy we adopt in this paper, is to assign a physical meaning to the purely 
formally defined normalisation factor and "densities" by looking at simple weighted walk 
problems for which we can compute the partition functions. The weights of the configura- 
tions depend on parameters which correspond to the rates of the stochastic processes and 
the partition function coincides with the normalisation factor defined in ((7J) if a common 
factor to all the cofactors is removed. In the next sections we illustrate this approach with 
the help of an example. We first consider a combinatorial problem which is interesting on 
its own. This is the one-transit walk model. We will compute its partition function and 
obtain the phase diagram of the model from the properties of the densities. We also show 
that the same partition function coincides with the normalisation factor of the TASEP 
model. 

3 The OTW model versus the TASEP 

The totally asymmetric simple exclusion process (TASEP) has grown to be one of the 
main theoretical models of nonequilibrium statistical physics. This is not only due to its 
simplicity and general applicablity, but also because its stationary probability distribution 
(SPDF) and other properties can be calculated exactly ^3EHE1E3E3- I n this paper 
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we show that this SPDF can be regarded as an equilibrium probability distribution of a 
simple model of a walk near an interface. Since the walk model is an equilibrium system, it 
can be described thermodynamically using standard methods. The phase behaviour of the 
TASEP can be explained in terms of adsorption transitions of the walk on the interface. 

In Section fci.il we introduce a model of a walk in the vicinity of a fixed interface. The 
walk is allowed to penetrate the interface once. Both ends of the walk are fixed but the 
point of penetration is free. An excess interface fugacity l/z\ is associated for contacts 
above the interface, and a fugacity l/zy for contacts below. For this model we are able 
to calculate the phase diagram exactly. The walk model without penetration was used as 
a simple model for polymer adsorption in [23]. After discussing the thermodynamics of 
the walk model, we show in Section f3.3l that it is closely related to the TASEP with open 
boundaries. More precisely, the statistical partition function Z{z\ 1 z 2 ) of the walk model 
is equal to the reduced normalisation of the stationary state of the TASEP if the interface 
fugacities in the walk model are equal to the in- and output rates of the two reservoirs. 

We show that the thermodynamic TASEP density p and current J are related to the 
contact densities p\ and p2, conjugate to z\ and Z2 respectively, as 

2/3 - 1 = 1 ~ P2 1 ~ Pi (n , 

J Z 2 Zl ' 

where 

1 duj 
lo = - lim -logZ n (zi,z 2 ), p i = -z i — . (12) 

rwoo n azi 

These equations show that one may derive the thermodynamic behaviour of quantities 
for a nonequilibrium model from those of an equilibrium model. We hope to be able to 
make contact between our approach and the formulation of a free energy functional for 
the TASEP from large deviation functions as adopted in [23GSI- We also would like to 
point out that the walk model has an appealing analogy with a continuous model for 
the dynamics of shocks in terms of which the TASEP phase diagram can be explained 
quantitatively [57] . 

3.1 The one-transit model 

Consider a statistical model of a restricted solid-on-solid (RSOS) path, also called Dyck 
path, on the rotated square lattice. Paths start at (0, 0) and end at (2n, 0), can only move 
in the North-East (NE) or in the South-East (SE) direction and cross the x-axis exactly 
once, see Figure ^ We associate energies — e± and —£2 to the returns (or contact points) 
of the path above and below the x-axis respectively. To make contact with Section |2J we 
implement this in the following way. A fugacity z\ = e £l l kT is given to each down step, 
and Z2 = e £2 / kT to each up step, except those ending on the x-axis. This model is directly 
related to the canonical model of 1281 . 
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-1 1 2n+l 

Figure 1: An example of an RSOS path starting at (0,0) and ending at (2n, 0) crossing 
the x-axis only once. 

By reflecting the last part of the RSOS path in the x-axis, it can be easily seen that 
the total number of possible paths is equal to the number of Dyck paths of length 2n. It 
is known, see e.g. |29| that the number of Dyck paths with p returns is given by Ballot 
numbers, 

pf2n — p—l\ p{2n — p — 1)\ 
n\ n — 1 J n\{n — p)\ 

The total number C n of Dyck paths of length 2n can be obtained by summing over p in 
(|13|) . or by noting that it is equal to the number of Dyck paths of length 2n + 2 with 
exactly one return, 

C n = B n p = B n+ n = — — ( ) , (14) 
p=i x 7 

which is the Catalan number. The partition function of the one-transit model is simply 
given by 

Z n (z 1 ,z 2 ) = (z 1 z 2 ) n Z n (z ll z 2 ), (15) 

where 

n p 

Z n ( Zl ,Z 2 ) = Y, Bn,P £ Z i qz 2 P+q - (16) 
p=0 g=0 

This can also be written in the following way, 

n 

Z n (z 1 ,z 2 ) = (z 1 z 2 ) n '^2z p (z 1 ,oo)Z n _ p (oo,z 2 ). (17) 

p=0 

This formula shows that we can also interpret our model as the combination of two contact 
models with a movable but impenetrable wall in between them at a random position, each 
position being equally probable. Equation 1)17(1 thus defines the partition function of an 
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annealed system, i.e. where the partition sum is averaged over the random position of the 
wall. 

The partition sum l|17[) is equal to the reduced normalisation of the totally asymmetric 
exclusion simple process (TASEP) [Toll 16U1 711 18U19| if the fugacities z\ and z 2 are replaced 
by its boundary in- and output rates. This result is important since it will allow us to 
associate the "densities" defined formally from the reduced normalization of the TASEP 
with the physical densities of the OTW model. In Section \3.'A\ we will go deeper into the 
relation between the OTW model and the TASEP, but before that we first describe the 
phase diagram of the OTW model. 



3.2 The phase diagram of the OTW model 

We define the grand potential per site for the gas of contacts as 



u> = — lim —\ogZ n . 

n— >oo n 



(18) 



Note that to get rid of spurious factors of 2 in subsequent formulas we divide by n instead 
of the system size is 2n. The potential uj can be easily calculated from (|17|) once we know 
the asymptotic properties of Z n (z, oo) = Z n (oo,z). This is well known, see e.g. [T51I25] . 
and can for example be derived from the differential equation it satisfies, 

(2n — IV 

- (1 - z)(l - 2z)Z' n (z, oo) + z(z + n(l - 2z) 2 )Z n (z, oo) = 2z 2 ± -=-A (19) 



n 



Analyzing the large n behaviour of this equation for the regions z > 1/2, z = 1/2 and 
z < 1/2 we immediately obtain 



Z n (z, OO) 



4" 



(1 - 2zf ^nV 2 
4™ 



l-2z 1 
1 - z z n (l - z) n 



z > 1/2 
z = 1/2 
z < 1/2. 



(20) 



The grand potential uj is given by minimizing over the position of the domain wall, and is 
therefore given by 

(21) 



uj(z 1 ,z 2 ) = -log 4z 1 z 2 + inf uj x (z 1 ,z 2 ), 

0<x<l 



where 



uj x (z 1 ,z 2 ) 



Zl ,z 2 >l/2 
x log (1 — z\) + (1 — x) log 4z2(l — z 2 ) elsewhere 



(22) 
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Prom (|21j) one finds the grand potential in all regions of the phase diagram, 



uj(zi,z 2 ) 



log4ziz 2 2l,Z2>l/2 

log z 2 + log(l - zi) z\ < 1/2, z 2 > zi 

log z\ + log(l - z 2 ) z 2 < 1/2, zi > z 2 



(23) 



We now turn to the calculation of the contact densities, which are the order parameters 
of the OTW model. From the definition of the walk model it immediately follows that 
the probabilities (aj) n and (6j) n to have a contact at site 2i above or below the x-axis, are 
given by, 



(bi)n 



{z\z 2 ) 
{z\z 2 ) 



Zi(zi,oo)Z n _i(zi,z 2 ) 

Z n {zi, z 2 ) 
^Z^zx, z 2 )Z n -i(co, z 2 ) 



Z n (z 1 ,z 2 ) 

We define the average number of contacts at x by 
and find in the thermodynamic limit n — > oo, 



(24) 
(25) 

(26) 



({a x ), (b x )) 



(0,0) 

(P(zi),0) 

(0, P (z 2 )) 



zi,z 2 > 1/2 

Z\ < 1/2, z 2 > z\ 

z 2 < 1/2, z\ > z 2 



(27) 



with 



P(z) 



1 - 2z 
1-z' 



(28) 



In this limit, these numbers are independent of x except on the line z\ = z 2 = z where we 
find 

((a x ),(b x )) = (p(z)(l-x),p(z)x). (29) 

The total number of contacts above and below are denoted by (a) and (6) respectively, 
and the corresponding thermodynamic densities can be calculated through derivatives of 
the grand potential, 



(a) „ duj 
a= hm — = 1 + Z\- — , 

n-^oc n OZ\ 



duj 
l + z 2 — . 

oz 2 



(30) 
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Note that lo is not everywhere differentiable. If uo is not differentiable in a point z* we 
define 

duo 1 / , Qlu , . Quj , \ 

2:* — (2*) = lim - (z* - e) — ( z * ~ e ) + ( z * + £ )^r\ z * + £ ) ) ■ ( 31 ) 



<9z 



+ o 2 



dz 



dz 



With this definition, (|30|) is valid everywhere. Because we average over the position of 
the domain wall, the densities a and b are not independent. Their values can be easily 
calculated and are given by, 



(a, b) 



(0,0) 

(P(*i),0) 

(0,pfo)) 
(p(z)/2,p(z)/2) 



zi,z 2 > 1/2 
zi < 1/2, z 2 > zi, 
z 2 < 1/2, zi > z 2 
z\= z<i = z < 1/2. 



(32) 



Notice that either both densities are equal or only one of the two does not vanish. This 
means that effectively one sees only one density. We thus find that for zi,z 2 > 1/2 the 
walk is entirely desorbed from the interface. When z\ < 1/2 and z 2 > z\ the walk is 
adsorbed above the interface, the contact density a is nonzero, while it is desorbed below 
and vice versa when z 2 < 1/2 and z\ > z 2 , see Figured 




desorbed 




desorbed 
desorbed 



adsorbed 



1/2 zi 

Figure 2: Phase diagram of the walk model. 
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The grand potential ()23|) is non- analytic at the lines z\ = 1/2 when z<i > z\ and 
Z2 = 1/2 when z\ > Z2- There is also a singularity at the line z\ = %i = z when z < 1/2. 
These lines therefore indicate phase boundaries. There is a first-order phase transition 
along the line z\ = Z2 = z for z < 1/2 along which the mirror symmetry of the system is 
spontaneously broken. The sum of the two densities r = a + b varies continuously across 
the line but their difference d = a — b is discontinuous. 

Above the line z<i = 1/2 (z\ > z%) the densities of contact points a = (a)/n and 
b = {b)/n vanish as n — > oo. Approaching the line z<i = 1/2 from above, the number of 
contacts (b) diverges like 

d) ~ ^n, (33) 

where we have used (|20j) . Using the same equation, on the critical line Z2 = 1/2 we get 

b ~ n" 1/2 = n-*, (34) 

and below the critical line 

b ~ 1 - 2z 2 = (1 - 2z 2 ) 1 /^ 1 . (35) 

A similar behaviour is obtained for the line z\ = 1/2 when z 2 > z\ and if we replace 
the density b by a. The critical behavior (|33|) and (|35|) as well as the finite-size scaling 
behavior (|34|) characterize a special surface phase transition [21] with a single critical 
exponent <p = 1/2. A similar exponent is found in other equilibrium problems with long 
range interactions [23] ■ The interest in discussing the phase diagram comes from the fact 
that it gives another physical interpretation of the phase transitions observed in TASEP. 

We conclude this section with a description of the model using the canonical ensemble. 
We now consider a and b as free parameters. The canonical free energy per site for given 
values of a and b can be calculated from the grand potential u)(zi, Z2), 

f(a,b) = sup((l - a)logzi + (1 - b)logz 2 +u(zi,z 2 )), (36) 

21,22 

from which we find 

/(a, b) = max{g(a,b),g(b,a)}, (37) 
g(a,b) = (l-a-6)log(l-a)-(2-a-6)log(2-o). (38) 

This result can be conveniently rewritten using r = a + b and d = a — b as 

/(a, b) = (1 - r) log f 1 - r -±^) - (2 - r) log ( 2 - . (39) 
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3.3 Connection with the totally asymmetric simple exclusion process 

In ^B] the stationary state of the TASEP was constructed using equivalent representations 
of the DEHP algebra. The reduced normalisation calculated using this method is equal to 
the partition sum of the OTW model, given by (|16|) . if the boundary rates of the TASEP 
are replaced by the contact fugacitites of the OTW model. In this section we show more 
precisely how the OTW model of Section 13.11 is related to the TASEP. 

Following an observation by Brak and Essam [2B] (see also see also [SO]) that the 
different equivalent representations of the DEHP algebra can be interpreted as transfer 
matrices for various lattice walk models, we construct a new representation which will give 
the transfer matrix for the OTW model. 

3.3.1 Transfer matrix formalism 

The OTW model can be described using a transfer matrix formalism. We will show that 
the partition function Z n {z\,z>i) can be written in the following form 



where T is the transfer matrix. We introduce a two-step transfer matrix T = T°T e , where 



The matrices D\ and E\ will act as transfer matrices for the walk above the x-axis, and 
D2 and E2 for the walk below the x-axis. The upper triangular form of T° ensures then 
that the walk can cross the x-axis only once. We will now describe the transfer matrices 
in detail. 

The matrix element (-Di)m for j > 2 is the weight of an edge from a point with height 
y = 2i — 2 to a point with height 2j — 3. The first column of D\ is auxiliary whose meaning 
will become clear later. Similarly, (-£>2)i? for j > 2 denotes the weight of an edge from a 
point with height y = 2 — 2i to a point with height 3 — 2j and the first column is again 
auxiliary. If the ket |n) represents the height n, the matrices D\, D2 and S are given in 
terms of projectors as, 



Z = (L\T n \R) 



(40) 




(41) 



CO 




(42) 



n=0 



00 




(43) 



n=0 



s 



z 2 |o)H + |o)(-i 



(44) 
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where \u) denotes an auxiliary ket vector. Explicitly, the matrices D\ and D2 are given 
by, 



/0 1 

110 

D x = 11 

1 

\0 

and the matrix S is given by, 



S 



0\ 



•7 



Do 



(xi 




Vo 



100 \ 
1100 
0110 
o o 1 '■ 
000 '••/ 



(45) 



(x 2 


1 








o\ 














































\o 











'••/ 



(46) 



The parameters x\ and X2 are arbitrary and will not enter the partition sum. For simplicity 
we could therefore set them to zero, but we will need them later for another reason (see 
equation (|57|l). 

The matrix element {E\)ij for i > 2 is the weight of an edge from a point with height 
y = 2i — 3 to a point with height 2j — 2. The first row of E\ is auxiliary. Similarly, {E-i)%j 
for i > 2 denotes the weight of an edge from a point with height y = 3 — 2i to a point 
with height 2 — 2j and its first row is again auxiliary. In terms of projectors, the matrices 
E\ and E2 are given by, 



71=1 



E x = x 3 |n)(0| +zf 1 |l)(0| + J^(|2n-l) + |2n + l))<2 

n=l 

00 

E 2 = z 2 - 1 |-l)(0|+J](|- 2n - 1 ) + |- 2n + 1 ))(- 2n l- 
Explicitly, they are given by, 

£1 



n 



(47) 
(48) 



/ x 3 











o\ 




( 











o\ 


H X 


1 















1 














1 


1 








E 2 = 





1 


1 














1 


1 













1 


1 





V 










'••/ 




V 













(49) 



Also here, the parameter X3 is arbitrary and will not enter the partition sum. 
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To indicate that walks can only start and end at height we furthermore define the 
vectors (L\ = (i(L|, 2 (£|) and \R) = (\R) 1 ,\R) 2 ), such that 

i(L| = (0| = (1,0,0,...), 2(^| = 0, 

(50) 

= |0> = (1,0,0,...), \R) 2 = |0> = (1,0,0,...). 

It is straightforward to check that the partition sum (|16|) of all walks of length 2n can be 
expressed as (|4*U|) . 

We end this section by defining the even and odd identity matrices for future conve- 
nience, 

'I?' e 



where 



L 2 



I°' e = ( {, jo,e ) , (51) 



1° = |0)(u| +^|2n)(2n- 1|, 1£ = |0)(u| + ^ |-2n)(-2n + 1|, (52) 

n=l n=l 

oo 00 

If = |«)(0| + ^|2n-l)<2n|, /f = |u)(0| + J2 \~ 2n + l )(~ 2n \- (53) 

n=l n=l 

3.3.2 The TASEP revisited 

The asymmetric simple exclusion process in continuous time is a particle hopping model 
with excluded volume in one dimension, where particles hop from the left to the right with 
rate 1. In the presence of open boundaries, the input rate of particles on the left of the 
system is a and the output rate on the right is (3, see Fig. |3 

a f3 



Figure 3: Sample TASEP configuration. Particles enter the system from the left with rate 
a and leave from the right with rate (3. Particles hop in the bulk from left to right with 
rate 1. 

If the Ti £ {0, 1} denotes the presence or absence of a particle, one would for ex- 
ample like to know the probability P(t\, t 2 , ■ ■ ■ , T n } to find a system in configuration 
{ T i, T 2, ■ ■ ■ ,T n } in the long time limit. In this limit, all these probabilities are stationary, 
and in ^H] this stationary state was calculated exactly. In JH] it was shown that this 
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solution can be conveniently expressed in a matrix product form, 

1 n 

P(ri,..., r n ) = ^{W\ H(nD + (1 - n)E)\V), (54) 
An i=i 

where the normalisation Z n is given by 

Z n = (W\(D + E) n \V), (55) 

and the matrices D and E, and the vectors (W\ and \V) are a representation of the 
so-called DEHP algebra, 

DE = D + E 

d\v) = hy) (56) 

(W\E = -(W\. 
a 

From the result ()16|) one may already have inferred that the partition sum of the 
OTW model is equal to the normalisation of the stationary state of the asymmetric sim- 
ple exclusion process (ASEP) with open boundaries ^H]- The fugacities Z\ and z 2 are 
then identified with the in- and output rates a and (3 respectively. Indeed, if we set the 
parameters x% to the values, 

xi = z^ 1 , x 2 = z 2 x , x 3 = zf 1 , (57) 

we find that the transfer matrices T° and T e , and the vectors (L\ and \R) constitute the 
following representation of the DEHP algebra (|56j) , 

D = T°I e , E = I°T e , \V) = \R), (W\ = (L\. (58) 

Various representations of |16j for the DEHP matrices were used as transfer matrices 
in |2H] to find bijections between several different path problems. Here we remark that 
this interpretation of the DEHP matrices allows to express a stationary nonequilibrium 
probability distribution in terms of an equilibrium distribution. Among other things, this 
has the consequence that the thermodynamics of the nonequilibrium model is prescribed 
by standard equilibrium thermodynamics. 

The TASEP can be formulated using a transition matrix, see e.g. |16j . The normalisa- 
tion as defined by ()55j) is not equal to that calculated from the cofactors of this transition 
matrix using the results of Section |^1 In this approach, for each system size, all the cofac- 
tors have a common factor which is a nontrivial polynomial. As hinted at in Section |2J we 
believe that this common factor will not give rise to additional singularities for positive 
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real rates in the thermodynamic limit. For example, upon introduction of inhomogeneities 
in the transition rates, the cofactors will no longer have a common factor. However, if 
these inhomogeneities are small enough the physical properties of the system should re- 
main the same. Moreover, we believe that a common factor is nonzero in the space of 
complex rates if the real parts of all the rates are positive. This is called the half-plane 
property, see e.g. We have checked this for small system sizes in the case of the 

TASEP. In the spirit of the Lee- Yang theory it implies that the common factor does not 
develop singularities in the thermodynamic limit for positive real rates and hence it will 
not influence the phase diagram, except perhaps at the origin. 

3.4 OTW-TASEP relation 

In this section we show how the TASEP current and density can be related to the equi- 
librium densities of the OTW model. 

3.4.1 Current 

The TASEP current operator is given by, 

j = (r°/ e )(7°r e ). (59) 

The average value, 

Jn,i = (Ji)n = ^(Lir-ijT"-*- 1 !^, (60) 

has the following meaning in the path problem. The two identity matrices in (|59|) have, 
above the x-axis, the effect of forcing an upstep between column 2i — 1 and 2i and a 
downstep between 2i and 2i + 1. Below the x axis they have the effect of forcing a 
downstep between column 2i — 1 and 2i and an upstep between 2i and 2i + 1. Therefore, 
Ji is the average number of paths that have a local maximum above or a local minimum 
below the x-axis between columns 2i — 1 and 2i + 1. The pieces of the path before and 
after these local extrema can be concatenated to obtain a path of length 2n — 2. Since 
the local extrema may occur at any height, we thus obtain all paths of length 2n — 2 and 
therefore 

Jn,i — Jn — 7> ' (^-Q 

independent of i. In the OTW model, the current corresponds to the pressure, since it is 
essentially the volume derivative of the grand potential. The value of the current 

J = lim J n , (62) 

n— »oo 
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in the various parts of the phase diagram is 

1/4 z 1 ,z 2 >l/2 

zi(l - z%) z\ < 1/2, z 2 > z\ 

z 2 (l - z 2 ) z 2 < 1/2, zi > z 2 



J 



(63) 



3.4.2 Density 

The contact operators can be given in terms of projectors, 

a-i = |1>2»— i<0|2i, k = |— 1>2»— 1<0| 2 < 
so that the contact number operators can be rewritten as 



a = ^fii, b = y^bj 



(64) 



(65) 



i=l 



i=\ 



The TASEP density operator fj also has an expression in terms of the projectors of the 
OTW model. The operator tj is obtained by putting the matrix I e instead of T e at position 
2i 



n = I e (2i). 



(66) 



This has the effect that between columns 2i — 1 and 2i each walk above the x-axis must 
go up. Walks below the x-axis must go down between these columns at all heights except 
y = —1, where it also may go up. From the result of Derrida et al. ^B]) or from a 
combinatorial argument it follows that the expectation value (ft) can be written as 



(n)n 



1 



Z n (zi,z 2 ) 



n—i—l 



z 2) 



p=0 



(67) 



where we have used (|14[) and H16|) . Using the expression for the expectation values of the 
contacts (oi) n and (bi) n , see eqs. (PH and l|25jl. we find 



n—i—l p j 

(Ti) n = ^2 IT ^ n ~i ^ Jn-l{bi-l)n-l, 

p=0 j=0 Z2 

and with the particle-hole symmetry of the TASEP this is equivalent to 

i-2 V 



p=0 3=0 



j Jn-l(o>i-l)n-l- 



(69) 
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Equations (|68(1 and ()69|) give the relations between the local densities of the OTW model 
and the TASEP. Combining (JHEJ) and (|MJ) we find 



n-i-l p , 1 1 s 

1+ C p JJ Jn-j + </n-l ( — (&i-l)n-l (Oi-l)n-l ] 

p=t_l 7=0 ^ 2 21 ' 



(70) 



In the bulk, the second term in the right hand side of (|7Ujl vanishes in the thermodynamic 
limit. We thus find that in each part of the phase diagram the following relation between 
the TASEP bulk density p and current J, and the equilibrium densities a and b is satsified, 



2p- 1 
J 



b a 

Z2 Zi 



(71) 



where z\ = a, %i = j3 and a and b are given by (f52|) . 



3.5 The partially asymmetric exclusion process 

The TASEP can be extended with a nonzero rate q for back-hopping. The resulting model 
is called the partially asymmetric simple exclusion process (PASEP). Exact results for the 
symmetric case (q = 1) are given in |32| and the stationary state of the general PASEP 
can be found using a matrix method as well |2UU21| . To the rate q will now correspond 
a "number of particles" N(q), not present in the TASEP. For the forward bias regime 
(q < 1) the phase structure is similar to the TASEP model, we are interested in a new 
phenomenon which occurs in the vicinity of q = 1. One finds 

N(q) = ^_ n + 0(logn), q < 1, a,/3>(l-g)/2 

1 " 1 (72) 

N(q) = ^n 2 +0(n), q>\ 

This implies a change of the volume when the transition rate q changes. We notice that 
the density p(q) = N(q)/n defined for q < 1 diverges for a finite value of q, namely q=l. 
For q > 1 we have to redefine the density as p(q) = N(q)/n 2 . This density turns out to 
be independent of the fugacity. A change of the volume was observed also in TASEP at 
the phase transition between the disordered state and the maximum current state. The 
latter phase transition could be interpreted as a special surface phase transition known in 
polymer physics (the number of "particles" were either proportional with the size of the 
system or independendent of the size of the system). Equation (|72|) describes a different 
phase transition since the number of "particles" are either proportional with the size of 
the system or with the square of the size of the system. We expect therefore that a 
simple extension of the OTW model could explain what one observes. This is indeed the 
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case. In the new model, a walk gets height dependent step weights. One can formulate 
a thermodynamical theory, analogous to that of a heterogeneous gas in a gravitational 
field [33] and obtain a partition function equal to the normalisation factor of For 
q > 1 the OTW model is genuinely two-dimensional, and has therefore a volume of order 
n 2 . For q < 1 the system undergoes a bulk phase transition and the only contributions to 
the grand potential now come from the surface, i.e. the system becomes one-dimensional. 
As in the TASEP, for q < 1 the system may undergo further phase transitions through 
enhancement of the surface chemical potentials. This is indeed the what happens and we 
find the adsorption-desorption transitions discussed in Section 13.21 

4 Conclusion 

In a previous paper ^1] we have shown that a properly chosen normalisation factor of the 
probability distribution function describing the stationary state of the "raise and peel" one- 
dimensional model is given by the partition function of the two-dimensional ice model with 
domain-wall boundary conditions — an equilibrium problem with nonlocal interactions. 
This connection was proven for small systems but there are good reasons |34j to believe 
that this conjecture is valid for any size of the system. It turns out that the same way 
to choose the normalisation factor was suggested in a general framework by Blythe and 
Evans |Hj and shown to be useful in order to use the Lee- Yang approach to non-equilibrium 
problems. This brought us to have a closer look at the problem. We have first noticed 
that the "number of particles" associated with various rates are non-decreasing functions of 
the rates seen as fugacities. This allows, as in equilibrium problems, to determine directly 
the phase diagram of a model, once the normalisation factor is known. This observation 
suggests, obviously, that one can try approximative approaches like finite-size scaling or 
power expansions to determine the nature of the phase transition in the case when the 
normalisation factor is not known exactly for all sizes. 

We have also shown, in the example of TASEP that, as in the "raise and peel" model, 
the normalisation factor can be understood as a partition function of a two-dimensional 
equilibrium model: the one-transit walk model. We also think that the correspondence 
between normalisation factors of one-dimensional stationary states and two-dimensional 
equilibrium problems with nonlocal interactions is of a more general validity. 
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A The normalisation as a homogeneous polynomial 



The normalisation of a stationary state of any stochastic (Markov) process can always be 
interpreted as a polynomial in the rates with positive coefficients, i.e. it has the form of 
a generating function. By the Cauchy-Schwartz inequality, the negative logarithm of this 
generating function is therefore convex and its derivatives with respect to the rates are 
proper "particle" numbers, i.e. the second derivatives are positive. 

The statement above is implied by the matrix-tree theorem |llU13llP2] . Here we show 
a simple proof which can be found in a slighly different version in |22| . 

Lemma 1 Let H be a matrix with off-diagonal elements H a i, = — r j and such that all 
columns add up to zero, ^2 a H a b = 0. Assume that H has a unique largest eigenvalue equal 
to 0. 

a) The cofactors X(a,b) are constant for each column, i.e. they do not depend on a. 

b) The eigenvector corresponding to the largest eigenvalue is a polynomial in the rates 
r a fc with positive coefficients. 

Proof. 

Let H(a, b) be the matrix corresponding to H with the oth row and 6th column re- 
moved. The cofactor X(a, b) is then defined by, 

X(a, b) = (-l) a+b det H(a, b). (73) 

Using row operations that do not change the determinant, and because of the special 
properties of H, it is not difficult to transform H(a,b) into H(b,b), and hence X(a,b) = 
X(b,b) for all a. 

Because the eigenvalue is unique and, 

= det 5 = Y J H ab X(a,b) = H ab X(b,b), (74) 

b b 

the elements of the eigenvector corresponding to the eigenvalue are given by the cofactors 
X(b, b). Each such cofactor is of the form 

n 

x(b,b)=j2 N ^--- 

,7Tf,_l,7r6 + l, . . . ,7T n J 

Il^oo (75) 

IT c=l 

where the sum is over any permutation ir = {tti, . . . , i^b-i-, ^b+i-, ■ ■ ■ ■> ^n} of {1, ... , b—1, b+ 
1, . . . , n} and N(ir) £ 7L. We now show that in fact N(tt) E {0, 1}, hence proving assertion 
b) of Lemma n 
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Let r Pa ^ a = 1 for a particular permutation p, and all other r ac = 0. From (|75jl we 
then see that N(p\, . . . , Pb-x, Pb+ii ■ ■ ■ > Pn) is the determinant of a matrix which we will 
call H (6, 6, p). If p a 7^ 6 for all a = 1, . . . , b — 1, b + 1, . . . , n, the columns of H (b, 6, p) all 
add up to zero and det H(b, b, p) = 0. If on the other hand p a = b for a particular a = a* , 
then H (b, b, p) contains zeros in the column corresponding to a* except for the diagonal 
element which is 1. By deleting the column and row of H(b,b,p) corresponding to a* we 
find again a matrix of the form of H (b, b, p) but with one dimension less. The result thus 
follows by expanding the determinant with respect to the column corresponding to a* and 
induction on n. 
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